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ABSTRACT 


This  report  presents  the  results  of  an  experimental  investigation  on  the  deformation  and  strength 
properties  of  several  particulate  cohesive  systems  under  quasi-static  compression  for  a  variety  of 
loading  histories  and  strain  rates  conditions.  These  systems  are:  (i)  silica  (silicon  dioxide)  powder, 
spherical  particles  having  a  mean  particle  size  of  63.5  pm;  (ii)  polyethylene  pellets,  cylindrical  in 
shape  (4  millimeters  in  height  by  3  millimeters  in  diameter),  and  (iii)  microscrystalline  cellulose, 
PH-105,  of  angular  shape  and  mean  particle  size  20pm.  The  experimental  activities  were 
complemented  with  model  development  in  the  framework  of  plasticity  and  viscoelasticity  theories. 
Comparison  between  theory  and  data  show  that  the  models  developed  describe  very  well  the  main 
features  of  the  behavior  observed  in  the  range  of  pressures  and  strain  rates  considered  in  the 
experimental  investigation.  A  new  model  for  calculating  the  resistance  to  penetration  into 
compressible  porous  viscoplastic  media  was  proposed.  The  influence  of  the  characteristics  of  the 
penetrator/target  interface,  impact  velocity,  target  mechanical  properties  and  nose  geometry  on  the 
resistance  to  penetration  were  investigated. 
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TASK  1:  EXPERIMENTAL  CHARACTERIZATION  OF  PARTICULATE 

SYSTEMS 


1.1  Introduction 

This  task  focus  is  on  the  experimental  characterization  of  the  compaction, 
deformation/flow  behavior,  and  strength  of  various  cohesive  particulate  systems.  A 
variety  of  index  tests  of  flow  ability  exist  such  as  the  Hosokawa  Powder  Tester,  or  a 
Copley  Tap  Density  Volumeter,  both  of  which  can  be  utilized  to  detennine  the  Hausner 
ratio  of  aerated  to  tapped  density,  an  often-used  index  of  flowability  (see  for  example 
Abdullah  and  Geldart,  1999).  Such  index  tests  are  often  useful  in  indicating  whether  one 
powder  will  flow  more  easily  than  another;  however,  they  provide  little  information  that 
can  be  used  in  a  comprehensive  model  of  the  flow  behavior.  For  many  industrial 
applications  a  Jenike  tester  has  become  the  default  standard  shear  testing  technique  used 
to  characterize  the  quasi-static  strength  of  powders  under  a  range  of  confining  loads.  This 
technique  suffers  from  several  shortcomings  including  sensitivity  to  the  filling  procedure 
and  induced  anisotropy  in  the  sample  through  the  initial  compaction  phase.  Furthennore, 
the  shear  is  non-unifonn  and  the  strain  deformation  is  quite  limited.  The  results  change 
by  as  much  as  a  factor  of  two  if  the  shearing  direction  is  changed  from  the  direction  of 
shear  during  the  compaction  phase  (Feise,  1998).  The  resulting  flow  function  is  not  a  true 
failure  surface  that  can  be  utilized  in  a  mathematical  constitutive  model  of  bulk 
deformation. 

A  simple,  yet  accurate  test  is  the  triaxial  compression  test  (Bishop  and  Henkel, 
1962).  In  the  triaxial  compression  test,  measurements  of  both  axial  and  volumetric  strains 
undergone  by  the  material  could  be  done  thus  allowing  the  characterization  of  stress- 
strain  response  under  axisymmetric  3-D  conditions.  Furthermore,  all  material  properties 
determined  by  direct  shear  testing  (ultimate  strength,  cohesion,  angle  of  internal  friction) 
could  also  be  estimated  from  triaxial  test  data,  with  comparable  results  (Kandala  and 
Puri,  2000). 

The  present  task  extends  the  current  knowledge  by  characterizing  the  quasistatic  behavior 
under  compression  for  a  variety  of  loading  histories  and  strain  rate  conditions.  The  focus 
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was  on  characterizing  the  behavior  under  dry  conditions  and  in  the  low  pressure  regime. 
To  gain  a  better  understanding  of  the  role  of  particle  size  distribution,  initial  porosity  and 
loading  history  on  compaction,  triaxial  compression  data  were  obtained  on  selected 
cohesive  systems:  (i)  silica  (silicon  dioxide,  floated  powder),  spherical  particles  having  a 
mean  particle  size  of  63.5  pm;  (ii)  polyethylene  pellets,  cylindrical  in  shape  (4 
millimeters  in  height  by  3  millimeters  in  diameter),  and  (iii)  microscrystalline  cellulose, 
PH-105,  of  angular  shape  and  mean  particle  size  20pm. 

Both  monotonic  and  cyclic  triaxial  compression  data  in  the  low  confining  pressure 
regime  are  reported.  By  performing  creep  tests  and  confined  compression  tests  at 
different  strain  rates,  the  time  influence  on  the  overall  behavior  as  well  the  strain  rate 
sensitivity  of  the  material  was  characterized. 

The  experimental  activities  presented  here  are  complemented  with  model  development 
for  the  systems  studied  in  Section  II. 

1.2  Experimental  Set-up 

Soil  mechanics  researchers  have  long  recognized  the  value  of  triaxial  testing  in  predicting 
the  stress-strain  and  failure  behavior  of  these  geologic  materials.  However,  the  classic 
soil  mechanics  triaxial  compression  apparatus  allows  measurement  of  volume  changes 
only  for  fully  saturated  conditions  (see  Bishop  and  Henkel  (1962)).  Yet,  the  potential  for 
significant  volume  changes  upon  shearing  under  dry  conditions  can  be  great,  since 
particles  can  form  loose  networks  of  contacts  which  can  move  as  rigid  bodies  and  cause 
dilation  of  the  bulk  material,  or,  can  collapse  to  denser  configurations.  A  few  researchers 
are  developing  equipment  capable  of  measuring  the  biaxial  and  triaxial  response  of 
particulate  systems  under  dry  conditions  ( e.g .  Abdel-Hadi  et  al.  (2000);  Arthur  (1985); 
Kandala  and  Puri  (2000);  Verwijs  et.  al.(2002);  etc.). 

The  equipment  that  was  used  for  testing  is  the  enhanced  triaxial  tester  developed  by 
Abdel-Hadi  (see  Abdel-Hadi  et  al.  (2000)).  The  setup  consists  of  five  major  components: 
(1)  the  triaxial  cell  that  contains  the  specimen  that  is  isolated  by  a  membrane  from  the 
confining  fluid  and  a  piston  that  applies  the  axial  load,  (2)  low  pressure  lines,  (3) 
electronic  volume  change  device,  (4)  data  acquisition  system,  and  (5)  a  universal  testing 
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machine.  The  major  enhancement  is  the  volume  change  device,  which  is  fully  automated, 
requires  no  human  intervention,  and  is  ideally  suited  for  computer  data  acquisition.  The 
elimination  of  cumbersome  pressure  lines  yields  improved  sensitivity,  accuracy  and 
reliability.  It  has  a  fast  response  time,  which  enables  it  to  detect  minute  or  abrupt  changes 
in  volume  during  testing.  A  schematic  diagram  of  the  testing  device  is  shown  in  Figure 
1.1. 


Triaxial  Cell 


Low  pressure 
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Air  pressure 

Electronic 
Volume  Change 
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/feedthrough 

Electrical  wiring 
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Figure  1.1:  Schematic  diagram  of  the  triaxial  tester  and  new  volume  change  device 
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A  typical  test  on  any  system  was  conducted  as  follows: 


(i)  Sample  preparation 

A  latex  membrane  of  0.3  mm  in  thickness  was  stretched  around  a  cylindrical  split 
mold  and  kept  in  place  by  applying  vacuum.  The  mold  was  sealed  by  rubber  O-rings  to  a 
pedestal  at  the  bottom,  and  to  a  cap  at  the  top.  The  material  was  poured  into  the  mold  in 
successive  layers  and  pre-compacted  using  an  integrated  compaction  device  that  can  be 
attached  to  the  split  mold.  This  new  device  is  shown  in  Figure  1.2.  Each  layer  was  thus 
pre-compacted  by  a  specific  amount  by  adjusting  the  drop  height  and  the  drop  weight  of 
the  pre-compaction  device.  In  this  way,  it  was  ensured  that  the  anisotropy  induced  by  the 
filling  procedure  is  negligible  and  that  all  the  specimens  are  pre-compacted  to 
approximately  the  same  density.  The  assemblage  was  then  placed  in  a  triaxial 
compression  cell  having  a  maximum  compression  capacity  of  3400  kPa.  The  cell  was 
filled  with  water,  and  then  placed  on  the  load  cell  of  an  Instron  5566  testing  machine. 

(ii)  Testing  protocol 

Any  test  consisted  of  2  stages:  (1)  hydrostatic  and  (2)  deviatoric.  In  the  hydrostatic  part  of 
the  test,  the  specimen  was  subjected  to  an  all-around  pressure  that  was  gradually 
increased  up  to  a  predetennined  level,  say  cr3.  The  quantities  measured  during  the 

hydrostatic  part  of  the  test  were:  (1)  the  applied  pressure,  (2)  the  change  in  the  volume  of 
the  sample,  which  was  determined  by  the  change  in  the  volume  of  water  used  to 
compress  the  specimen  (Abdel-Hadi  et  al.  (2000)).  Once  the  desired  level  of  confining 
pressure  cr3  was  reached,  the  lateral  pressure  was  kept  constant  while  an  axial  load  was 
applied  to  the  sample  through  the  loading  piston  (deviatoric  part  of  the  test).  Both  axial 
and  volumetric  strains  were  continuously  recorded.  All  tests  were  conducted  in  a 
thermostatically  controlled  air-conditioned  laboratory  with  minimal  ambient  changes,  the 
recorded  temperatures  and  relative  humidity  being  25°C  and  45%,  respectively. 
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Figure  1.2:  Device  and  split  mold  for  controlled  pre-compaction  of  the  particulate 
systems  prior  testing. 
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1.3.  Experimental  Results  on  Silica  Powder 
1.3.1.  Material  Characteristics 

The  silica  powder  (silicon  dioxide,  floated  powder)  under  study  is  monosized  with  a 
mean  particle  size  of  63.5  pm  or  240  mesh.  The  geometrical  configuration  and  surface 
structure  of  the  system  is  shown  in  an  SEM  (scanning  electron  microscope)  micro-graph 
(see  Figure  1.3).  The  particle  density  was  about  2.3  g/cm  and  the  initial  density  of  the 
powder  was  about  1.3  g/cm  .  Although  this  system  is  fairly  well  characterized  in  terms  of 
surface  chemistry,  the  available  information  concerning  its  mechanical  behavior  is 
limited.  Another  rationale  for  studying  this  system  rather  than  sand  is  because  it  is 
monosized  and  all  particles  have  the  same  geometry  so  we  could  separate  the  effects  of 
particle  size  distribution  and  geometry  from  that  due  to  confinement  on  the  compaction 
of  the  material.  Here  and  throughout  the  text  compressive  stresses  and  strains  are 
considered  positive. 


Figure  1.3:  SEM  micrograph  of  monosized  silica  of  of  63.5  pm  mean  particle  size  (or 
240  mesh). 
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Monotonic  triaxial  compression  tests  under  lateral  confining  pressure  cr3  in  the  range  20 
kPa  to  350  kPa  and  cyclic  triaxial  compression  tests  consisting  of  several  loading-creep¬ 
unloading  and  reloading  cycles  under  the  same  confining  pressures  were  conducted. 
These  tests  allowed  the  measurement  of  both  the  deformability  and  strength 
characteristics.  Measurements  of  the  volume  changes  of  the  sample  as  a  function  of  the 
applied  stress  were  used  to  detennine  the  evolution  of  the  bulk  modulus  i.e.  compaction 
characteristics  of  the  material  as  well  as  the  yield  behavior  under  3-D  stress  conditions. 

1.3.2.  Monotonic  Test  Results  on  Silica  Powder 

To  gain  a  fundamental  understanding  of  the  compaction  properties  of  the  material,  the 
experimental  program  included  monotonic  triaxial  compression  tests  under  lateral 
confining  pressure,  cri in  the  range  20  kPa  to  350  kPa;  the  crosshead  displacement  rate 

was  0.1  mm/min.  A  slip-stick  phenomenon  is  observed  (sudden  reduction  in  material 
strength)  at  the  beginning  of  the  tests  for  the  very  low  confining  pressures. 

Typical  results  of  a  hydrostatic  compression  test  up  to  350  kPa  are  shown  in  Figure  1.4. 
As  an  example  of  a  deviatoric  test,  the  stress-strain  curves  (i.e.  axial  and  volumetric 
strain  vs.  the  stress  difference  cr,  -  cr3 )  obtained  under  a  confining  pressure  cr3  =  93  kPa 

are  shown  in  Figures  1.5  (a)  and  (b).  The  initial  density  of  the  sample  was  pi  =1.33 
g/cm  .  The  stress  versus  axial  strain  curve  presents  3  regions:  (1)  quasi-linear  elastic,  (2) 
concave  towards  the  strain  axis,  (3)  plateau.  In  the  first  part  (up  to  200  kPa),  the  behavior 
can  be  considered  to  be  elastic  (the  irreversible  strains  are  small  in  comparison  to  the 
elastic,  recoverable  strains).  In  the  second  region,  the  plastic  deformation  undergone  by 
the  material  is  significant  and  the  value  of  the  yield  limit  continuously  increases  up  to 
failure  at  <jd  =  <rl  -  cr3  =  320  kPa.  The  existence  of  a  plateau  shows  that  the  material  has 
viscous  (viscoplastic)  properties,  that  is,  it  deforms  even  under  constant  stress. 
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■  Experimental  data 
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Figure  1.4:  Pressure  vs.  volumetric  strain  response  under  hydrostatic  compression  for 
silica. 
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Figure  1.5(a):  Deviatoric  stress  vs.  axial  strain  at  93  kPa  of  confining  pressure. 
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Figure  1.5(b):  Deviatoric  stress  vs.  volumetric  strain  at  93  kPa  of  confining  pressure. 


The  stress-volumetric  strain  curve,  indicates  that  the  material  undergoes  further 
compaction  under  deviatoric  conditions,  and  then  dilates  (see  Figure  1.5(b)).  The 
existence  of  those  two  regimes  of  volumetric  behavior  was  observed  in  all  tests.  Figure 
1 .6  shows  the  volumetric  profiles  obtained  in  tests  under  confinement  levels  of  40  kPa, 
70  kPa,  and  93  kPa. 
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Figure  1.6:  Volumetric  strain  profiles  at  confining  pressures  of  93  kPa,  40  kPa,  and  70 
kPa 

It  is  clear  that  the  usual  decoupling  between  the  deviatoric  and  hydrostatic  responses  that 
is  usual  in  constitutive  models  belonging  to  the  commercial  finite  element  software  for 
transient  analyses  such  as  LS-DYNA  cannot  be  considered  in  the  case  of  silica. 
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Note  that  the  stress  value  for  which  the  curve  a 1  -  cr3  versus  the  volumetric 
deformation sv  changes  sign  depends  on  the  confining  pressure.  These  critical  points 

define  the  boundary  between  the  compressible  and  dilatant  regimes.  As  is  the  case  with 
most  particulate  systems,  the  strength  of  the  material  is  an  increasing  function  of  the 
applied  confining  pressure.  Figure  1.7  shows  the  deformed  specimen  at  the  end  of  a 
deviatoric  test  under  70  kPa  confining  pressure. 


Figure  1.7:  Triaxial  compression  cell  with  deformed  specimen. 
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1.3.3.  Cyclic  Test  Results  on  Silica  Powder 


To  characterize  the  influence  of  the  loading  history  and  time  effects  on  the  response  of 
silica,  cyclic  confined  compression  tests  were  also  perfonned  for  the  same  confining 
pressures  as  in  the  monotonic  tests.  The  results  obtained  in  a  test  under  27  kPa  confining 
pressure  are  shown  in  Figures  1.8  (a)  to  1.8  (c).  The  test  consisted  of  five  loading¬ 
unloading  and  reloading  cycles.  Before  passing  from  loading  to  unloading  the  load  was 
held  constant  for  about  20  minutes  in  order  to  separate  possible  viscous  effects  from 
unloading.  Then,  unloading  and  further  reloading  was  performed.  From  the  unloading 
slopes  of  a j  -  cr3  versus  ex,  the  Young  modulus  E  was  detennined.  A  dependence  of  E 

CT  _|_  2^ 

on  the  mean  stress  p  =  — 5 — — — -  is  observed.  From  the  unloading  slopes  of  cr,  -  cr3 

versus  sv  curves,  the  bulk  modulus  K  was  evaluated.  As  an  example,  in  Table  2  are 

given  the  values  of  the  elastic  moduli  detennined  from  data  obtained  in  a  cyclic  triaxial 
compression  test  under  a  confining  pressure  of  27  kPa. 


Cycle  Number 

E  (MPa) 

K  (MPa) 

1 

63.500 

168 

2 

57 

161 

3 

63 

181 

4 

70.4 

195.35 

5 

97.04 

254.34 

Table  1.1.  Elastic  constants  as  detennined  from  data  obtained  in  a  cyclic  triaxial 
compression  test  on  silica  under  confining  pressure  of  27  kPa. 
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(kPa) 


Note  also  a  change  in  the  creep  direction  at  the  transition  from  the  compaction  regime  to 


Axial  Strain  (mm/mm) 

the  dilatancy  regime.  This  effect  can  be  clearly  seen  in  a  zoom  of  the  stress-volumetric 
strain  curve  shown  in  Figure  1.8  (c). 

Figure  1.8(a):  Deviatoric  stress  vs.  axial  strain  curve  in  a  cyclic  test  under  27  kPa 
confining  pressure. 


14 


(kPa) 


120 


Figure  1.8(b):  Deviatoric  stress  vs.  volumetric  strain  in  a  cyclic  test  under  27  kPa 
confining  pressure. 
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Figure  1.8(c):  Zoom  in  of  the  deviatoric  stress  vs.  volumetric  strain  curve  in  a  cyclic  test 
under  27  kPa  confining  pressure  showing  the  change  in  creep  direction  during  transition 
between  the  compressibility  and  dilatancy  regime. 
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1.4.  Experimental  Results  on  Polyethylene  Pellets 


1.4.1.  Material  Characteristics 

In  this  subsection,  we  present  the  results  of  a  series  of  triaxial  compression  tests 
conducted  on  commercially  available  dry  PE  pellets.  The  composition  of  the  pellets  is: 
polyethylene  >  99%  and  talc  <  0. 1%- 1 .0%.  Their  appearance  is  translucent,  white  and 
soft.  In  general,  the  pellets  are  uniform  and  cylindrical  in  shape  (4  millimeters  in  height 
by  3  millimeters  in  diameter).  The  pellet  density  is  0.91  g/cnr ,  while  the  initial  density  of 

'j 

the  PE  system  varied  from  0.59-0.63  g/cm  . 

1.4.2.  Monotonic  Test  Results  on  Polyethylene  Pellets 

Typical  results  of  a  hydrostatic  compression  test  up  to  350  kPa  on  polyethylene  pellets 
(PE)  are  shown  in  Figure  1.9.  It  can  be  concluded  that  the  behavior  of  PE  pellets  under 
isotropic  compression  is  linearly  elastic  within  the  pressure  range  considered. 

A  series  of  monotonic  triaxial  compression  tests  were  perfonned  for  confining  pressure  in 
the  range  of  14.7  kPa  to  35  kPa;  the  crosshead  displacement  rate  was  of  0.1  mm/min.  For 
illustration  purposes,  we  present  the  test  results  of  the  experiments  under  confining 
pressures  at  the  extremes  of  the  pressure  range.  In  Figure  1.10  (a)  the  stress-axial  strain 
curves  are  shown.  It  can  be  noted  that  the  stress-strain  response  is  non-linear.  For  a  given 
confining  pressure,  the  material  underwent  elastic  deformation  followed  by  irreversible 
deformation  and  failure.  The  corresponding  stress-volumetric  strain  curves  (Figure  1.10 
(b))  show  two  distinct  regimes  of  behavior:  (1)  a  compressibility  regime,  followed  by  (2) 
a  dilation  regime.  The  occurrence  of  both  regimes  of  volumetric  behavior  was  observed 
in  all  tests  in  the  given  range  of  confining  pressures.  Figures  1.10(a)  and  (b)  also  show 
the  effect  of  the  confining  pressure  on  the  overall  response,  i.e.  the  higher  the  confining 
pressure,  the  higher  the  yield  limit,  dilatancy  threshold,  and  strength.  An  increase  of 
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Mean  Pressure  (kPa) 


confining  pressure  of  138%  from  14.7  to  35  kPa  resulted  in  an  increase  of  123%  in 
compressive  strength  and  of  120%  in  the  dilatancy  threshold  level.  However,  for  any 
given  confining  pressure  the  dilatancy  threshold  level  corresponds  to  about  75%  of  the 


Volumetric  Strain  (mm3/mm3) 

compressive  strength. 

Figure  1.9:  Pressure  vs.  volumetric  strain  response  under  hydrostatic  compression  for  PE 
pellets. 
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Figure  1.10(a):  Stress  -  axial  strain  curves  showing  the  effect  of  confining  pressure  on  the 
behavior  of  polyethylene  pellets. 


Princ 


Volumetric  Strain  (mm3/mm3) 


Figure  1.10(b):  Stress-volumetric  strain  curve  showing  compressibility  followed  by 
dilatancy. 


Strain  rate  sensitivity 

Figures  1.11  show  the  results  of  a  standard  triaxial  compression  test  at  cr3=  35  kPa  and 

0.1  and  10  mm/min  strain  rates,  respectively.  A  clear  rate  sensitivity  of  the  overall 
response  can  be  noted:  an  increase  in  strain  rate  results  in  steeper  slope  of  the  stress-strain 
curve,  higher  dilatancy  threshold,  more  dilatancy,  and  higher  strength  (a  1.25  fold 
increase  in  strength  as  a  result  of  a  hundred-fold  increase  in  the  strain  rate). 


Axial  Strain  (mm/mm) 

Figure  1.11(a):  Stress-axial  strain  curves  in  monotonic  triaxial  compression  test  at  cr3  = 
14  kPa  showing  the  strain  rate  sensitivity  of  polyethylene  pellets; 
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Volumetric  Strain  (mm3/mm3) 

Figure  1.1 1(b):  Stress-  volumetric  strain  curves  in  monotonic  triaxial  compression  test  at 
cr3  =  14  kPa  showing  the  strain  rate  sensitivity  of  polyethylene  pellets. 
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Figure  1.12:  Comparison  between  the  volumetric  stress-strain  behavior  of  coated  and 
non-coated  polyethylene  pellets  showing  that  the  coated  system  exhibits  much  more 
dilatancy  and  thus  improved  flowability. 


Coated  vs.  non-coated  pellets  behavior. 


The  stress-strain  response  of  a  coated  ethylene  and  styrene  copolymer  sample  in  pellet 
form  was  investigated  and  compared  to  the  behavior  of  non-coated  polyethylene  pellets 
(see  Figure  1.12).  The  results  show  that  although  the  dilatancy  threshold  is  the  same  for 
both  systems,  the  coated  sample  exhibits  much  more  dilatancy.  Thus,  the  enhanced 
triaxial  tester  provides  an  excellent  means  for  comparing  the  flowability  of  different 
systems. 
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1.4.3.  Cyclic  Test  Results  on  PE  pellets 


To  characterize  the  influence  of  the  loading  history  and  time  on  the  response  of  the 
material,  cyclic  confined  compression  tests  were  also  perfonned  for  the  same  confining 
pressures  as  in  the  monotonic  tests.  As  an  example,  in  Figures  1.13  (a)  and  (b)  are  shown 
the  results  of  a  test  corresponding  to  cr3  =  35  kPa.  The  test  consisted  of  four  loading- 

creep-  unloading  -reloading  cycles.  Note  that  before  unloading,  the  axial  load  was  held 
constant.  The  material  deformed  by  short-term  creep,  which  appeared  to  stabilize  in  15- 
20  minutes.  Then,  partial  unloading  and  further  reloading  was  perfonned.  A  change  in 
creep  direction  at  the  transition  from  the  compressive  to  the  dilatancy  regime  was 
observed(see  Figure  1.14).  Lastly,  note  the  absence  of  hysteresis  loops  and  that  the 
unloading  curves  are  quasi-linear.  The  Young’s  modulus  E  was  detennined  from  the 
slopes  of  (cr,  -  cr3)  versus  axial  strain  (s] )  curves  during  unloading.  A  dependence  of  E 

on  stress  state  was  observed  (see  Figure  1.13  (a)).  Similarly,  the  bulk  modulus,  K  was 
evaluated  from  the  unloading  slopes  of  deviatoric  stress  vs.  volumetric  strain  curves  and 
its  variation  with  the  stress  level  is  indicated  in  Figure  1.13  (b).  A  variation  of  both  elastic 
parameters  with  stress  was  observed  in  all  tests. The  mean  arithmetic  values  calculated 
from  all  tests  was:  E  =  27555  kPa  and  K=  17975  kPa,  respectively.  It  should  be  noted  that 
the  experimental  procedure  (short-creep  prior  to  unloading)  allowed  for  a  very  good 
separation  of  viscous  effects  from  unloading,  and  consequently  increased  accuracy  in  the 
estimate  of  the  elastic  properties  of  the  material.  More  details  concerning  this 
experimental  procedure  and  its  application  to  a  variety  of  geologic  materials  can  be  found 
in  Cristescu  and  Hunche  (1998). 
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Figure  1.13(a):  Stress-axial  strain  curves  in  cyclic  triaxial  compression  test  at  cr3=  35 

kPa;  before  each  partial  unloading,  axial  load  was  kept  constant;  Young  modulus  E 
values  computed  from  the  unloading  slopes  are  indicated. 
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Figure  1.13(b):  Stress-  volumetric  strain  curves  in  cyclic  triaxial  compression  test  at  cr3  = 
35  kPa;  before  each  partial  unloading,  axial  load  was  kept  constant;  computed  values  of 
the  bulk  modulus  K  are  indicated 
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Figure  1.14:  Detail  of  the  stress-volumetric  strain  curve  at  cr3=  35  kPa  showing  the 
change  in  creep  direction  as  the  material  transitions  from  compressible  to  dilatant 
behavior. 
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Stress-path  and  short-term  time  effects :  Based  on  the  cyclic  test  results,  it  can  be 
concluded  that  the  material  displays  elastic/viscoplastic  behavior.  The  question,  which 
arises,  is  whether  the  dominant  properties  of  PE  pellets  are  elastic/plastic  or 
elastic/viscoplastic.  In  other  words,  what  is  the  importance  of  short-term  time  effects  on 
the  response?  A  comparison  between  the  material  response  for  monotonic  and  cyclic 
loading  is  shown  in  Figure  1.15.  In  both  tests  the  crosshead  displacement  rate  during 
loading  was  of  0.1  mm/mln.  It  appears  that  short-tenn  creep  effects  can  be  neglected  and 
that  an  elastic/plastic  model  could  describe  accurately  the  main  features  of  the  behavior: 
non-linearity  of  the  overall  stress-strain  response,  hardening  volumetric 
compressibility/dilatancy,  and  an  influence  of  the  stress  invariants  on  the  behavior. 


Axial  Strain  (mm/mm) 

Figure  1.15:  Comparison  between  the  stress-axial  strain  curves  at  cr3=  35  kPa  for 
mono  tonic  and  cyclic  loading  paths. 
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1.4.4  Long-term  Creep  Test  on  PE  pellets 


A  long-term  creep  test  under  confined  conditions  (lateral  confining  pressure  of  14.7  kPa) 
was  also  conducted.  The  axial  load  was  kept  constant  at  160  N  (cTj=  52.53  kPa)  for  24 

hours  and  both  the  axial  and  volumetric  strains  were  continuously  monitored.  Figure  1.16 
(a)  shows  that  the  stabilization  of  the  axial  strain  took  place  in  about  20  minutes.  A  very 
interesting  evolution  in  time  of  the  volumetric  strain  was  recorded  (see  Figure  1.16  (b)). 
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Figure  1.16(a):  Evolution  of  the  axial  strain  with  time  in  a  creep  test  over  24  hours  at 
cr3  =  14  kPa  and  cr,  =  52.53  kPa  showing  that  stabilization  is  reached  within  20  minutes. 
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Figure  1.16(b):  Volumetric  strain  versus  time  curve  showing  that  volumetric  creep  does 
not  stabilize  within  the  duration  of  the  test  (24  hours  at  cr3  =  14  kPa  and  cr,  =  52.53  kPa 


Although,  it  appeared  that  the  volumetric  strain  stabilizes  in  20  minutes,  a  sharp  increase 
in  creep  rate  was  observed  between  l-10hrs,  followed  by  a  stabilization  plateau,  and 
again  a  change  in  creep  rate  after  about  20  hrs  from  the  beginning  of  the  creep  stage. 
Stabilization  was  not  reached  within  the  test  duration  (24  hrs).  The  test  results  show  a 
marked  difference  between  axial  and  volumetric  creep  behavior.  Furthermore,  the  volume 
viscosity  is  not  constant,  but  depends  on  the  current  density  (or  level  of  compaction). 
Clearly,  such  effects  cannot  be  simulated  in  the  framework  of  the  classic  theory  of 
viscoplasticity  based  on  the  concept  of  overstress. 
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1.5.  Experimental  Results  on  microscrystalline  cellulose 


1.5.1.  Material  Characteristics 

The  material  studied  is  microcrystalline  cellulose,  PH-105,  (C6Hi0O5)n,  mean  particle  size 
20pm.  It  has  a  bulk  density  of  about  0.25g/cm3  and  it  is  highly  compressible.  This 
microcrystalline  cellulose  (MCC)  is  monosized  and  has  particles  spherical  in  shape.  It 
was  tested  to  provide  information  on  the  role  of  particle  size  into  the  deformation  and 
strength  through  comparison  with  the  test  results  on  silica  (see  Section  1. 3.). 


1.5.2.  Test  Results  on  microcrystalline  cellulose 

Hydrostatic  cyclic  tests  with  several  loading-creep-unloading-partial  reloading  cycles  up 
to  a  pressure  of  275.8  kPa  were  perfonned.  The  average  value  of  the  bulk  modulus  was 
determined  to  be  of  245  kPa.  A  series  of  monotonic  and  triaxial  compression  tests  were 
conducted  for  confining  pressures  ranging  from  68.9  kPa  to  275.8  kPa.  As  an  example, 
in  Figure  1.17  are  shown  the  deviatoric  stress-volumetric  strain  curves  obtained  for 
confining  pressures  of  68.9,  137.8,  and  206.7  kPa,  respectively.  For  any  given  level  of 
confinement,  the  material  showed  only  compressibility,  the  final  volumetric  defonnation 
increasing  with  the  confinement  (up  to  4  times  when  lateral  confinement  is  increased 
from  68.9  to  206.7  kPa).  For  any  given  level  of  confinement,  a  monotonic  increasing 
variation  of  the  bulk  modulus  with  the  applied  stress  deviator  was  observed  (see  Figure 
1.18).  The  effects  of  strain  rate  of  compaction  were  also  investigated  for  crosshead  speed 
varying  from  0.1  mm/min  to  100  mm/min.  The  stress-volumetric  strain  curves  obtained 
under  the  same  confinement  but  at  an  increased  strain  rate  show  that  the  higher  the  strain 
rate,  the  higher  is  the  degree  of  compaction  (see  Figures  1.19). 
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Figure  1.17:  Deviatoric  stress  vs.  volumetric  strain  curves  for  microcrystalline  cellulose 
under  confining  pressures  of  69  kPa  ,138  kPa,  and  207  kPa. 
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Bulk  Moduli,  K  (kPa) 


Figure  1.18:  Variation  of  the  bulk  modulus  K  for  microcrystalline  cellulose  as  a  function 
of  deviatoric  stress  showing  increased  compaction  under  shear  conditions. 


Volumetric  Strain,  ev 


Figure  1.19(a):  Strain  rate  effects  on  compaction  as  observed  from  the  variation  of 
volume  in  a  deviatoric  test  for  microcrystalline  cellulose  at  a  confining  pressures  of  69 
kPa  for  crosshead  speed  of  lOOmm/min. 
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Figure  1.19(b):  Strain  rate  effects  on  compaction  as  observed  from  the  variation  of 
volume  in  a  deviatoric  test  for  microcrystalline  cellulose  at  a  confining  pressures  of  69 
kPa  for  crosshead  speed  of  lOmm/min 


References  Section  I 


Abdel-Hadi,  A.  2004.  Development  of  a  new  technique  for  measuring  volume  change  of 
dry  particulate  systems  under  very  low  confining  pressures,  Ph.D.  Dissertation,  Univ  of 
Florida,  Gainesville. 

Abdullah,  E.C.  and  Geldart,  D.,  1999.  The  use  of  bulk  density  measurements  as 
flowability  indicators.  Powder  Technol.  102,  151-165. 

Arthur,  J.R.F.,  Dunstan,  T.,  Enstad,  G.G.,  1985.  Determination  of  the  Flow  Function  by 
Means  of  a  Cubic  Plane  Strain  Tester.  Int.  J.  Bulk  Solids  Storage  Silos,  1  (2),  99-128. 
Bishop,  A.  W.,  Henkel,  D.  J.,  1962.  The  Measurements  of  Soil  Properties  in  Triaxal  Test, 
2nd  ed.,  Edward  Arnold,  London. 

Cristescu,  N.D.,  Hunsche,  U.,  1998.  Time  Effects  in  Rock  Mechanics.  John  Wiley  & 
Sons,  Ltd.,  Chichester. 

Feise,  H.J.,  1998.  A  Review  of  Induced  Anisotropy  and  Steady-State  Flow  in  Powders, 
Powder  Technol.  98,  191-200. 

Kandala,  R.N.,  Puri,  V.M.,  2000.  Measurement  of  cohesion  and  angle  of  internal  friction 
using  cubical  triaxial  tester  and  comparison  with  computer  controlled  shear  cell,  Particul. 
Sci.  Tehnol.  18(2),  71-88. 

Verwijs,  M.J.,  Janssen,  R.J.M.,  Scarlett,  B.,  2002.  Influence  of  the  intermediate  principal 
stress  during  powder  shear  in  a  biaxial  tester.  In:  Proceedings  of  the  World  Congress  of 
Technology  of  Particles.  Sydney,  Australia,  July  2002  (CD-ROM). 


36 


TASK  2:  CONSTITUTIVE  MODELING  OF  THE  RESPONSE  OF  THE  PARTICULATE 

SYSTEMS 


II.  1  Introduction 

This  task  focus  is  on  the  development  of  robust  constitutive  models  at  the  continuum  level  that  can 
describe  the  behavior  of  particulate  systems  investigated  over  a  wide  range  of  stress  and  strain 
conditions.  A  large  body  of  literature  exists  in  which  particulate  systems  are  modeled  as  a  plastic 
continuum  using  soil  mechanics  models.  Critical  state  models  such  as  the  Cam-clay  model  (e.g. 
Roscoe  and  Burland,  1968)  have  been  applied  to  cohesive  powders  (e.g.  Puri  et  al.  (1995)).  Double¬ 
surface  models  such  as  the  Drucker-Prager/cap  elasto-plastic  model  have  been  used  to  describe 
ceramic  powder  stress-strain  behavior  (see  Gethin  et  al.  (1994),  Aydin  et  al.  (1996);  Chtourou  et  al. 
(2002),  etc.). 

Continuum  level  models  formulated  in  the  framework  of  viscoplasticity  theory  have  been  used  for 
describing  short-term  time  effects  on  the  quasi-static  deformation  of  cohesive  systems  (e.g.  Cazacu 
et.  al.  (1997),  Jin  and  Cristescu  (1998),  Abdel-Hadi  et.  al.  (2002)  etc.).  These  models  are  capable  of 
simulating  nonlinear  elastic  response  (stress  dependency  of  the  elastic  parameters),  hardening, 
compressible/dilatant  volumetric  behavior,  creep  and  relaxation  at  room  temperature. 

All  the  models  discussed  so  far  neglect  the  effect  of  the  third  stress  invariant  on  the  behavior.  Yet, 
there  is  experimental  evidence  showing  that  all  three  principal  stresses,  ax ,  a2,  a3  (or  equivalently 

all  three  stress  invariants)  influence  the  deformation  and  failure  of  particulate  systems  (see  for 
example,  Bardett  and  Lode,  1990;  Peric  and  Ayari,  2000;  Alawaji  et  al.  (1992);  Lade  and  Kim 
(1995),  Saxena  et  al.  (1988),  Verwijs  et  al. (2002),  etc.).  Several  models  for  granular  materials  that 
make  use  of  yield  functions,  flow  potentials,  and  failure  conditions  involving  the  third  stress 
invariant  were  proposed  by  Lade  and  co-workers  (e.g.  Lade  (1975);  Lade  (1977);  Lade  and  Kim 
(1995),  etc).  Zhang  et  al.  (1986)  demonstrated  that  the  elastoplastic  model  developed  by  Lade 
(1977)  for  cohensionless  sands  adequately  describes  stress-strain  behavior  for  granular  agricultural 
materials. 

As  part  of  this  research  effort,  we  have  proposed  an  elasto-plastic  model  for  the  description 
of  the  inviscid  behavior  of  the  PE  system.  Concerning  the  silica  system,  we  examined  the 
possibility  of  modeling  the  inviscid  behavior  using  an  elasto-plastic  approach  and  the  time 
dependent  behavior  using  a  viscoelastic  type  approach.  For  both  systems  and  both  models,  a  general 
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procedure  for  the  determination  of  the  constitutive  parameters  based  on  a  minimum  set  of  triaxial 
compression  test  data  is  provided. 

In  section  II. 2,  we  present  the  elasto-plastic  model  adopted  for  PE  system.  It  is  shown  that 
the  inviscid  behavior  can  be  described  with  accuracy  using  a  single  surface,  Lade  (1977)  type 
elasto-plastic  theory.  A  robust  and  computationally  inexpensive  algorithm,  the  convex  cutting  plane 
algorithm  (Ortiz  and  Simo,  1986)  is  used  to  implement  the  elastic-plastic  model.  We  conclude  with 
a  comparison  between  the  model  prediction  and  data,  and  suggestions  for  future  work. 

Application  of  the  proposed  elastic/plastic  formulation  to  the  silica  system  is  presented  in 
section  II. 2.  Then,  we  proceed  with  the  presentation  of  the  general  form  of  the  proposed  constitutive 
viscoelastic  equation.  Finally,  through  comparison  between  model  predictions  and  data  we  show 
that  the  viscoelastic  model  provides  a  more  realistic  description  of  the  material's  response  than  the 
elastic/plastic  one. 

II.2  Elastic/Plastic  constitutive  model  for  PE  system 

From  the  test  results  conducted,  it  appears  that  short-term  creep  effects  can  be  neglected  and  that  an 
elastoplastic  model  could  describe  accurately  the  dominant  properties  of  the  PE  pellets.  To  reduce 
the  complexity  of  the  model,  we  assume  that  the  elastic  behavior  is  expressed  by  Hooke’s  law  with 
constant  moduli:  E  =  27555  kPa  and  K=  17975  kPa,  which  correspond  to  the  arithmetic  mean  of  the 
experimental  values  from  all  performed  tests.  For  the  sake  of  simplicity,  in  the  description  of  the 
plastic  behavior,  we  assume  an  isotropic  hardening  law,  i.e.  the  yield  surface  expands  uniformly 
and  that  the  degree  of  hardening  is  independent  on  the  stress  path.  Since  the  material  displays  both 
compressible  and  dilatant  behavior,  the  hardening  parameter  should  be  chosen  such  as  to  reflect  this 
pattern  of  the  plastic  volumetric  deformation.  Thus,  for  PE  pellets  we  assume  that  the  hardening 
parameter  is  the  plastic  work  per  unit  volume  defined  as: 


(1) 


where  ^  represents  the  Cauchy  stress  tensor  and  els’’  is  the  plastic  strain  increment.  The  plastic 
work  (1)  can  be  written  as 


Wr=\llde[ll+\&ltde'lf , 


(l.a) 
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where  /,  is  the  first  invariant  of  the  stress  tensor  ( +a2  +  cr3),  d.£y  is  the  plastic  volumetric 

strain  increment  and  “prime”  stands  for  deviator.  The  second  term  on  the  right-hand  side  of  (la)  is 
always  positive  and  represents  the  energy  input  needed  for  change  in  shape.  The  first  tenn  on  the 
right-hand  side  is  the  energy  related  to  plastic  volume  change  that  is  stored  during  compaction  and 
released  during  volumetric  expansion  of  the  system.  Because  Lade’s  elastic/plastic  theories  account 
for  the  dependence  of  the  plastic  deformation  and  strength  on  all  stress  invariants,  and  the 
hardening  parameter  is  the  plastic  work,  for  PE  pellets  we  use  the  expressions  of  the  failure,  yield 
function,  and  plastic  potential  proposed  in  Lade  (1977).  Thus,  the  failure  condition  for  PE  pellets  is 
assumed  to  be  of  the  fonn: 
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where  /3  =  cr,  <j2  cr3  stands  for  the  third  stress  invariant,  pa  is  the  atmospheric  pressure  measured  in 
the  same  units  as  the  stress,  while  m  and  ^  are  material  parameters.  The  numerical  values  of  these 
coefficients  are:  m  =  0.12  and  //,  =  21.4.  These  values  were  determined  using  strength  data  obtained 

in  standard  triaxial  compression  tests  (see  Section  1.3)  in  conjunction  with  (2).  In  the  principal 
stress  space,  the  shape  of  the  failure  surface  defined  by  (2)  is  conical,  with  the  apex  of  the  cone 
being  at  the  origin  of  the  stress  axes.  Since  the  hypothesis  of  isotropic  hardening  was  made,  under 
continuous  loading  the  yield  surface  expands  symmetrically  about  the  hydrostatic  axis,  its  ultimate 
position  being  the  failure  surface  (2).  Hence,  the  yield  function  is  given  by: 
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The  specific  expression  of  the  work-hardening  law  for  PE  pellets  can  be  determined  empirically  by 
calculating  the  experimental  values  of  the  plastic  work  in  a  standard  triaxial  compression  test  and 


plotting  its  variation  with  /  = 
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(see  Figure  2.1).  An  exponential  variation  as 


proposed  by  Lade  (1977)  was  found  to  approximate  well  PE  data,  i.e.: 
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where  the  coefficients  a,  b,  and  c  are  functions  of  the  confining  pressure  (see  Figure  2.1  showing 
the  relationship  between  f  and  the  plastic  work  at  14.7  and  35  kPa  confining  pressure, 

respectively). 


Figure  2.1:  Variation  of  the  plastic  work  normalized  by  the  atmospheric  pressure  with  the  stress 
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at  <ji  =  14.7  and  35  kPa,  respectively 
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For  PE  pellets,  a  linear  variation  of  c  with  the  confining  pressure  was  observed: 


c  =  k  +  d;  ■ 


a, 


P, 


(5) 


where  k  =  14047  and  g  =  -0.0447.  Since  the  ultimate  position  of  the  yield  surface  is  the  failure 
surface,  for  a  given  value  of  the  confining  pressure,  fp  ( Wpeak  j  =rj1 ,  where  Wpeak  is  the  value  of  the 
plastic  work  at  the  peak  deviatoric  stress.  Thus,  using  (4)  it  follows  that: 


a  =  P  i 
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where  e  is  the  basis  of  the  natural  logarithm.  The  variation  of  Wpeak  with  the  confining  pressure  can 


be  approximated  by: 


wr‘=n-p.- 
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(7) 


where  /  =  0.125  and  n  =  0.977.  Next,  by  making  use  of  (6)  and  (7),  the  laws  of  variation  of  the 
coefficients  a  and  b  with  confining  pressure  can  be  obtained. 

The  formulation  of  the  model  is  completed  by  providing  the  plastic  potential.  If  the  plastic  potential 
coincides  with  the  yield  function,  then 


del  / def  =(dfp/do-1)/(dfp/dcr3),  (8) 

dsp  and  dsp  being  the  plastic  axial  strain  and  plastic  lateral  strain  increments,  respectively,  cr,  the 
axial  applied  stress,  and  cr3  the  given  confining  pressure.  However,  for  PE  pellets  as  for  most 

particulate  systems,  (8)  is  not  in  agreement  with  the  experimental  data,  thus  the  flow  rule  ought  to 
be  nonassociated.  Following  Lade  (1977),  we  consider  for  the  plastic  potential  a  mathematical  form 
similar  to  that  of  the  yield  function,  i.e.: 
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where  r/2  is  a  constant  for  given  values  of  f  and  cr3 .  Note  that  if  r/2  were  equal  to  f  at  all  stress 
levels,  the  plastic  potential  would  coincide  with  the  yield  function  (2).  Thus, 
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where  dX  is  the  plastic  factor.  Next,  using  the  flow  rule  (10),  r/2  can  be  expressed  as: 
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where  vp  = - .  For  any  given  value  of  the  confining  pressure,  //2 ,  computed  using  (11),  was 

ds{ 

found  to  have  a  linear  variation  with  /  ;  the  slope  could  be  considered  constant,  whereas  the 

intercept  depends  on  the  confining  pressure.  A  simple  expression  of  the  type  proposed  by  Lade 
(1977)  was  found  to  approximate  well  PE  data: 


Th=s-fp+r 


(12) 


where  s  =  0.033,  r  =  -3.265,  and  t  =  -1.554.  Again,  since  the  experimental  values  of  rj2  are  smaller 
than  /  ,  /  does  not  coincide  with  gp  i.e.  the  flow  rule  is  nonassociated.  Note  that  negative  values 
of  rj2  correspond  to  compressive  volumetric  behavior  whereas  positive  values  correspond  to 
expansive  plastic  behavior  (Lade,  1977). 
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II.3.  Numerical  Implementation  and  Comparison  with  experimental  results 


The  increased  complexity  of  the  constitutive  models,  justified  by  the  need  to  reproduce  with 
improved  accuracy  material  behavior  under  general  loading,  deformation,  and  strain  rate  conditions, 
result  in  increased  complexity  of  the  numerical  algorithms.  In  an  effort  to  obtain  a  better  global 
accuracy  and  stability  of  the  numerical  solution,  algorithms  that  can  be  applied  to  a  specific  class  of 
models  are  being  developed  (see  for  example,  Ahadi  and  Krenk  (2003)  for  a  3D  generalization  of 
Cam-Clay  models  or  Borja  et  al.  (2003)  for  smooth  three  invariant  representations  of  Mohr- 
Coulomb  model). 

In  this  work,  a  robust  and  computationally  inexpensive  method,  the  convex  cutting  plane  (CCP) 
algorithm,  proposed  by  Ortiz  and  Simo  (1986)  was  used  to  integrate  the  single-surface  elastoplastic 
model  presented  in  the  previous  section.  It  is  an  implicit  procedure,  which  is  computationally 
attractive  because  it  can  be  linearized  in  closed  form.  The  nonnality  condition  is  enforced  at  the 
initial  (known)  iteration  and  the  algorithm  does  not  require  the  gradient  of  the  flow  rule  or 
hardening  law.  The  flow  chart  for  the  procedure  applied  to  the  elastoplastic  Lade  type  model  is 
given  below. 


1 .  Initialize  the  iteration  counter  step  k=0. 

2.  Initialize  the  plastic  strain  £PJ  )  :=  £p  ,  the  plastic  work  Wp  i<0)  =  Wpn  and  plastic  factor 

:=  0 . 

3.  Compute  C,  the  matrix  of  elastic  moduli,  possibly  depending  on  the  initial  stress 
<jn  through  its  invariants 

4.  While  k<k  max  repeat  steps  5-14. 

5.  Compute  the  trial  elastic  stress  <Ti+1(i’  =  C  (sn+1  —  SPn+[  ]. 

6.  Compute  [Ii (k>  1^  h<k>]  the  invariants  of  the  trial  elastic  stress  tensor  a n , ,  lk\ 


7.  Compute  the  increase  in  hardening  parameter  fp  „  +/k>  as:  f 


p  n+l 


~b-W 

{k)=a-e  pn+l 


Wpn+1 
V  Pa  J 


8.  Compute  the  plastic  moduli,  D  n+/k>  as 

1  1 


Dn+\k)  =a-e~bWpn+l 


-b- 


W, 


3 


p  n+l 


K  Pa  ) 


1 

+  -' 
c 


r  n 

[Paj 


■  w. 


p  n+l 


9.  Compute  the  yield  function  / 


(*) 


=Aa- 


+i 


» Wn 


+1 


r  ]3  w 
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-  W 
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10.  Is  fn+\k)  <  TOLx  true? 

The  condition  verities  that  cr,I+1<A)  is  on  or  under  the  yield  surface  of  hardening  parameter 

fp(n+i/k)  and  TOLi  is  a  tolerance,  accounting  for  rounding  errors. 

11.  If  YES,  the  final  new  stress  is  equal  to  the  elastic  trial  stress.  END  AND  EXIT. 

12.  If  NOT,  plastic  correction  is  needed,  as  follows: 


13. Obtain  increments  to  plastic  factor:  A2/L(| 


/„ 


(*) 


+i 


( (*)V 


Vn- 


da 


C 


¥jk'  a/M  n  m 

Scr  df 


n+ 1 


14.Update  plastic  strain,  plastic  factor,  and  state  variable: 


(*+ 1) 


c  —  cp  <A)  -i-  A2;W  — ;,+1 

£n+\  —  £n+ 1  +  ^  Aj+1  - 


3fn- 


do 


a-) 


(i-) 


r-  (t+l)  _  y’ 
p  n+1  J  i 


(k)  ,  a2;Wd 
^„+i  +A  /tn+1Dn+1 


15.  Update  iteration  counter  k=k+l  and  return  to  step  4. 


Figures  2.2  (a)  and  2.2  (b)  show  comparisons  between  numerical  simulations  and  data  obtained  in 
monotonic  tests  at  14.7  and  35  kPa  confining  pressures.  Note  that  the  model  describes  well  the 
overall  stress-strain  behavior  observed  in  the  tests.  The  main  features  of  the  volumetric  response, 
particularly  the  increase  in  the  value  of  the  dilatancy  threshold  level  with  increasing  confining 
pressure  are  very  well  simulated.  The  very  good  agreement  between  the  experimental  pressure- 
volumetric  strain  curve  under  isotropic  compression  and  the  linear  elastic  response  curve 
corresponding  to  K=  17975  kPa  (see  Figure  2.3)  shows  that  this  average  value  of  the  bulk  modulus 
is  an  accurate  measure  of  the  compressibility  of  the  material. 
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6 


Principal  Stress  Difference  (kPa) 


Figure  2.2  (b):  Comparison  between  theoretical  and  experimental  stress-volumetric  strain  curves  at 
cr3  =  14  kPa  and  35  kPa,  respectively; 
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Figure  2.3:  Comparison  between  theoretical  and  experimental  stress  strain  curves  for  hydrostatic 
compression  on  PE  pellets 


II.4.  Constitutive  Modeling  of  the  Behavior  of  the  Silica  System 
II.4.1  Elastic/plastic  Approach 


Data  shows  that  even  in  the  elastic  regime  the  behavior  is  strongly  nonlinear  (see  Section  1.3.2  and 
1.3.3).  Thus,  the  elastic  behavior  should  be  modeled  using  a  differential  Hooke’s  law  with  variable 
elastic  moduli.  To  reproduce  the  stress  dependency  of  the  Young  modulus  E,  we  propose  a  law  of 
variation  with  the  mean  pressure  of  the  form: 


-<5- 

E  =  Ex  +  y  ■  pa-e 


(13) 


where  /,  denotes  the  first  stress  invariant  (7,  =  <j1  +  cr2  +  cr3 )  ,  Ex (=  2.5  •  105 kPa)  is  an 
asymptotic  value  corresponding  to  very  high  pressures,  pa  is  the  atmospheric  pressure,  while  y  and 
8  are  material  constants.  For  silica,  we  found  that  y  =  -2.0017  •  105  and  8  =  8.3  •  10~4 . 
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Similarly,  the  variation  of  the  bulk  modulus  K  with  the  mean  stress  was  approximated  by: 


-P— 

Pa 


K  =  Kx+if/- pa-e  p°  (14) 

where  Kx  =2.5-10 5  kPa  ^  =  -1.399-105  and  /?  =  1.  21-10  3.  The  comparison  between  the 
theoretical  variation  of  the  bulk  modulus  with  the  mean  pressure  and  data  is  shown  in  Figure  2.4. 


Figure  2.4:  Comparison  between  the  experimental  and  theoretical  variation  of  the  bulk  modulus 

with  the  mean  pressure  for  silica  system. 

First,  we  interpreted  the  data  in  terms  of  the  Mohr-Coulomb  criterion.  In  this  model  it  is  assumed 
that  the  material  is  perfectly  plastic,  i.e.  the  yield  locus  is  fixed  in  the  stress  space.  The  yield  surface 
coincides  with  the  failure  surface: 

r/=c  +  cr„tan^  (15) 

where  z ,  and  crn  are  the  shear  and  normal  stresses  on  the  failure  plane,  respectively;  c  is  the 
cohesion,  and  (p  is  the  angle  of  internal  friction.  We  estimated  the  cohesion  and  the  angle  of  internal 
friction  from  compression  data.  We  found  c  =  0.56kPa  and  cp  =  33° .  Figure  2.5  shows  a 
comparison  between  model  and  data.  The  Drucker  Prager  criterion  was  also  used  to  interpret  the 
data.  The  advantages  of  using  Drucker-Prager  over  Mohr-  Coulomb  are:  (l)it  is  expressed  in  terms 
of  stress  invariants  and  thus  can  be  applied  to  describe  failure  surface  under  any  loading  conditions, 
(2)  it  has  been  proven  to  accurately  describe  the  strength  characteristics  of  granular  media.  The 
expression  of  Drucker-Prager  criterion  is: 

^h~D  =  a-Il+k  (16) 
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where  I\  is  the  first  stress  invariant,  /2D  the  second  invariant  of  the  stress  deviator,  a  and  k  are 
material  parameters.  For  silica  we  found  a  =  0.3197  and  k  =  1.304  . 


Figure  2.5:  Mohr-Coulomb  Yield  Surface  for  silica  system. 


We  have  also  modeled  the  observed  behavior  using  the  elastic/plastic  model  developed  by  Lade  and 
co-workers.  In  this  theory  the  total  strain  increment  ds  due  to  a  stress  increment  da  is  divided  in 
three  components:  an  elastic  strain  component  dse ,  a  plastic  collapse  strain  increment  component 
dsc ,  and  a  plastic  expansive  strain  increment  component  dsp ,  such  that: 

deM=del+del+e’  (17) 


The  elastic  strain  is  computed  using  a  non-linear  Hooke’s  law.  The  plastic  collapse  strain  is  given 
by  a  plastic  stress-strain  theory  involving  a  cap-type  yield  surface  while  the  plastic  expansive 
strains  are  computed  using  a  stress-strain  theory,  which  involves  a  conical  yield  surface  with  apex 
at  the  origin  of  the  stress  space.  The  collapse  behavior,  which  occurs  under  hydrostatic  conditions  is 
modeled  by  the  yield  function: 


,  r  w  y  2 

fc  =  7 -  pa 

\k  P a  ) 

where  W  is  the  irreversible  stress  work  and  fc 


(18) 

/|2  +  2  •  I22 ,  with  7j  and  /2  being  the  first  and 


second  invariant  of  the  stress  tensor.  For  our  system,  k  =  0.0122  and  d  =  0.8  . 


The  expansive  behavior  is  governed  by  the  following  yield  function: 

fP  = 


(1?  „  ^ 

(l) 

-4 — 27 

yh  7 

[Pa, 

(19) 
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while  the  failure  condition  is  given  by:  /  =  77, .  For  silica,  we  found  thatm  =  0.088  and^t  =  1.82  . 

Hydrostatic  compression  of  an  isotropic  powder  system  results  in  equal  linear  strains  in  the  three 
principal  directions.  For  this  condition,  the  direction  of  the  strain  increments  concides  with  the 
hydrostatic  axis.  Thus,  the  plastic  potential  must  be  equal  to  the  spherical  yield  cap  (fc ).  In  Figure 
2.6.  (a)  is  shown  the  location  of  the  yield  cap  relative  to  the  conical  failure  surface  for  silica.  The 
plastic  potential  for  expansive  strains  is  given  by: 


Sp  = 


f 

(„  ^ 

m  ^ 

27  +  r]  2 

Pa 

V 

l  A  ) 

) 

■h 


(20) 


For  given  values  of  f  and  cr3 ,  r]2  is  a  constant.  Negative  values  of  r]2  correspond  to  plastic 

volumetric  strains,  which  are  compressive,  and  positive  values  of  i]2  correspond  to  plastic  expansive 
strains.  In  Figure  2.6b  the  yield  surface  of  silica  in  the  principal  stress  space  is  also  shown. 


Figure  2.6:  (a)  Calculated  Failure  Surface  (-x-)  and  yield  cap  in  triaxial  plane. (b)  trace  of  the  yield 
surface  in  octahedral  plane  for  silica  system. 
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II.4.2  Numerical  Implementation 


Hydrostatic  Response 

The  Closest  Point  Projection  and  the  Cutting-Plane  algorithms  were  both  implemented  in  MatLab. 
A  parameter  sensitivity  study  was  done  on  the  constants  involved  in  the  expression  of  the  bulk 
modulus  and  hardening  law.  It  appears  that  the  closest  agreement  between  simulation  and 
experimental  results  is  obtained  when  the  Cutting-Plane  algorithm  is  used,  an  exponential  type  of 
law  of  variation  for  K  and  a  hardening  constant  d  =  1  (see  Figure  2.7). 


Deviatoric  Response 

Comparison  between  experimental  and  simulation  results  for  silica  show  that  the  model  describes 
particularly  well  the  effect  of  confining  pressure  on  the  volumetric  strains.  The  non-linearity,  the 
stress-path  dependency  and  the  gradual  change  from  compressible  to  dilatant  behavior  were 
correctly  modeled.  It  was  thus  demonstrated  that  the  theory  is  applicable  to  silica  for  general  3-D 
conditions.  As  an  example  in  Figures  2.8  and  2.9  is  shown  a  comparison  between  experimental  data 
and  simulation  results  using  the  model  for  a  confining  pressure  of  70.6  kPa.  Although  the  model 
overestimates  the  data,  the  transition  from  compressibility  to  dilatancy  is  well  reproduced. 
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Figure  2.7:  Mean  stress  vs.  volumetric  strain  for  silica  system.  Effects  of  hardening  constant  d. 
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- SmUation  (70.6  \<Pa) 

■  Experiments  (70.6  kRa) 


Axial  Strain  (mmfrrm) 


Figure  2.8:  Simulated  and  experimental  deviatoric  stress  vs.  axial  strain  at  70.6  kPa  for  silica 
powder. 
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- Simulation  (70.6  kPa) 

A  Experiments  (70.6  kPa) 


Figure  2.9:  Simulated  and  experimental  deviatoric  stress  vs.  volumetric  strain  curve  corresponding 
to  a  confining  pressure  of  70.6  kPa  confinement  for  silica  powder. 


II.4.3  Viscoelastic  Modeling 

In  this  section,  we  present  a  novel  quasi-linear  viscoelastic  model  aiming  at  describing  besides  time 
effects  such  as  creep,  stress  relaxation,  and  strain  rate  effects,  the  fundamental  mechanical  property 
exhibited  by  particulate  systems  that  is  to  be  either  compressible  or  dilatant  depending  on  the  stress 
state. 
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The  general  form  of  the  proposed  constitutive  equation  is: 


T) 

ev=-kv(ev-H(p,q ))  +  — 

K 

s'  =  -k(s'  -G(p,q)cr')  + 

1  Lr 


(21) 

(22) 


where  s'  denote  the  deviator  of  the  strain  tensor  and  sv  stand  for  the  volumetric  strain.  We  assume 
that  the  instantaneous  response  is  elastic  and  is  characterized  by  K  and  G,  the  dynamic  bulk 
modulus  and  shear  modulus,  respectively.  In  equations  (21)-(22),  H(p,q)  and  G(p,q)  define  the 
stabilization  boundaries  for  the  volumetric  and  shearing  strains  and  are  assumed  to  depend  on  stress 


tr<j 


through  p  =  -^~  and  q  =  .  —tr(c r')2,  the  first  invariant  of  the  stress  tensor  cr  and  the  second 


invariant  of  the  stress  deviator  cr' .  The  volumetric  and  shearing  viscosities,  denoted  by  k  and  kv , 

are  considered  to  be  constant.  The  main  features  of  our  model  that  distinguishes  it  from  the  classical 
viscoelastic  ones  are:  (i)  the  volumetric  part  of  the  strain  depends  on  both  the  deviatoric  (shear) 
stress  components  and  mean  pressure,  (ii)  the  deviatoric  part  of  the  strain  is  also  influenced  by  the 
volumetric  defonnability  characteristics.  We  will  show  that  by  a  judicious  choice  of  the 
mathematical  expressions  of  the  relaxation  boundaries,  all  the  observed  characteristics  of  the 
deformation  of  particulate  systems  could  be  described  with  high  accuracy. 


Application  of  the  constitutive  model  to  silica 

For  a  given  particulate  system,  the  specific  mathematical  expressions  of  the  constitutive  functions 
involved  in  the  proposed  model  could  be  determined  from  the  quasi-static  confined  triaxial 
compression  data.  Indeed,  the  monotonic  response  curves  at  a  very  low  strain  rate  could  be 
considered  to  be  close  to  the  stabilization  boundaries.  Under  confined  triaxial  compression 
conditions,  silica  exhibits  both  compressibility  and  dilatancy  and  the  ai  -  cr3  vs.  sv  curve  show  a 

plateau  for  stresses  close  to  failure  (see  Section  I,  Figure  1.5(b)).  To  reproduce  these  characteristics, 
we  approximate  the  relaxation  boundary  for  the  volumetric  strain  with  a  hyperbola  that  passes 
through  the  origin 
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ad 


H(p,q)  =  —  — — 

c  c[q-a) 

(23) 

and  of  asymptotes:  q  =  a  and  q  =  csv  +  d  ,  where  a,  c,  and  d  are  functions  of  the  stress  state.  For  a 
given  confining  pressure  cr3  =  constant,  a  represents  the  horizontal  asymptote  toward  which  tends 
the  cr,  -  cr3  vs.  sv  curve  whereas  c  is  the  double  of  the  slope  of  the  tangent  at  the  origin  of  the  same 
curve.  Furthermore, 


d  =  (q„-4  .tan;«  +  l  (24) 

a  tan“  a  - 1 

where  c  =  tan  2a  and  qd  denotes  the  dilatancy  threshold.  For  silica,  we  found  that  the  variation  of 
the  functions  a,  c,  and  qd  with  the  stress  state  is  very  well  described  by 


a  =  a(cr3)  =  a{<J3  +  a2 

(25.1) 

C  —  —  C|(T3  +  CjCJ 3  +  c3 

(25.2) 

< Id  =dd{^)=  Viai  +  V2^3  +  ^3 

(25.3) 

where  a,=  3.42,  a2=  68.16  kPa,  c,=  -0.1785  kPa'1,  c2=  38.836,  c3=  -707.82  kPa, 
v,=  -0.0228  kPa'1,  v2=  5.67,  and  v3=  -85.6  kPa.  Figure  2.10  shows  a  comparison  between  model 
predictions  and  data  for  the  confining  pressure  cr3  =  93  kPa.  According  to  the  constitutive  model 
under  confined  conditions,  the  evolution  equation  for  shear  strains  is: 

£i-£3=qg(p>q)  (26) 

where  ex  denotes  the  axial  strain,  £\  is  the  radial  strain,  p  =  (  at  +  2cr3 )  /  3  and  q  =  cr,  -  <r3 .  Thus, 
the  specific  expression  of  g(p,q)  can  be  obtained  by  approximating  the  stabilization  boundary. 
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For  silica,  it  is  found  that  under  a  given  confining  pressure  cr3  =constant,  the  following  polynomial 
law  fits  the  data  well 


g(q)  =  mlq5  +m2  (27) 

where  mx  and  m2  are  constants.  The  variation  of  mx  and  m2  with  the  confining  pressure  is  well 
approximated  by 


/«!  =  mx  (0-3 )  =  — j-  (28.1) 

t2  +  CT3 

m2  =  .v,  exp  (.S' 2  cr3 )  (28.2) 

where  tx  =  3.133  10'6,  t2  =  8.1 107  105,  ^=0.0022  and  ^,=  -0.0255. 


Figure  2.10:  Comparison  between  model  prediction  and  data  under  cr3  =  93  kPa 
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Figure  2.11:  Comparison  between  the  theoretical  variation  of  the  parameter  mx  and  data  (symbols). 


II.5.  Summary  and  conclusions 

An  enhanced  triaxial  testing  device  has  been  used  to  investigate  the  quasi-static  behavior  of  several 
cohesive  systems  under  dry  conditions:  (i)  silica  (silicon  dioxide,  floated  powder)  of  spherical 
shaped  particles  having  a  mean  particle  size  of  63.5  pm;  (ii)  polyethylene  pellets,  cylindrical  in 
shape  (4  millimeters  in  height  by  3  millimeters  in  diameter),  and  (iii)  microscrystalline  cellulose, 
PH-105,  of  angular  shape  and  mean  particle  size  20pm. 

Based  on  the  test  results,  we  have  adopted  an  elastic/plastic  approach  to  characterize  the  inviscid 
behavior.  The  model  describes  particularly  well  the  effect  of  confining  pressure  on  the  volumetric 
strains.  The  non-linearity,  stress-path  dependency  and  the  gradual  change  from  compressible  to 
dilatant  behavior  were  correctly  modeled.  It  was  thus  demonstrated  that  the  theory  is  applicable  to 
silica  for  general  3-D  conditions. 

To  describe  time  effects  on  the  material  behavior,  a  general  3  -D  quasi-linear  viscoelastic 
constitutive  model  was  proposed.  Contrary  to  existing  rate -type  constitutive  equations,  this  model 
predicts  a  dependency  of  the  volumetric  strain  on  both  the  deviatoric  (shear)  stress  components  and 
mean  pressure.  Furthermore,  it  was  shown  that  by  a  judicious  choice  of  the  mathematical 
expressions  of  the  relaxation  boundaries  both  compressibility  and  dilatancy  as  well  as  instantaneous 
failure  could  be  well  predicted.  A  procedure  for  the  determination  of  the  material  parameters  based 
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on  a  minimal  set  of  data  was  provided.  A  very  good  comparison  was  found  between  the  model 
prediction  and  data. 

Based  on  the  test  results  on  the  PE  pellets  it  can  be  concluded  that  under  hydrostatic  compression 
the  stress-strain  response  is  linear  and  elastic.  Under  deviatoric  conditions,  the  main  characteristics 
of  the  response  are:  non-linearity,  stress-path  dependency,  and  gradual  change  from  compressibility 
to  dilatancy  with  increasing  deviatoric  stresses.  The  time  influence  on  the  behavior  was  revealed  in 
creep  tests  as  well  as  in  tests  performed  at  different  strain  rates.  It  was  found  that  the  time 
dependent  properties  depend  on  the  loading  history.  However,  short-tenn  time  effects  can  be 
neglected  and  the  behavior  can  be  described  accurately  using  a  single  surface,  Lade  type  elastic- 
plastic  theory.  It  was  demonstrated  that  this  theory  captures  the  change  in  volumetric  deformation 
mechanisms  with  the  increase  in  the  applied  deviatoric  stress.  The  increase  in  the  value  of  the 
dilatancy  threshold  level  with  increasing  confining  pressure  is  particularly  well  described. 

Long-term  (over  24hrs)  creep  test  results  showed  a  marked  difference  between  axial  and  volumetric 
creep  behavior.  Furthennore,  the  volume  viscosity  is  not  constant,  but  depends  on  the  current 
density  (or  level  of  compaction).  Clearly,  such  effects  cannot  be  described  in  the  framework  of 
viscoplasticity  theory,  which  is  based  on  the  concept  of  overstress  and  where  the  yield  and  potential 
have  the  same  mathematical  expressions  as  the  yield  and  potential  for  inviscid  behavior.  However, 
a  viscoplastic  fluid  type  formulation  should  be  more  appropriate.  Additional  rheological 
measurements  would  be  required  for  the  development  and  calibration  of  such  a  model. 
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TASK  3:  ANALYSIS  OF  THE  STEADY-STATE  FLOW  OF  A  COMPRESSIBLE 

POROUS  VISCOPLASTIC  MEDIUM 


III.l.  Introduction 

The  focus  of  this  task  was  on  describing  steady  flow  of  a  compressible  porous  viscoplastic 
medium  over  a  penetrator.  The  influence  of  the  characteristics  of  the  penetrator/target 
interface,  impact  velocity,  target  mechanical  properties  and  nose  geometry  on  the  resistance 
to  penetration  is  investigated. 

Penetration  mechanics  has  a  long  and  rich  history.  A  comprehensive  review  of  empirical 
equations  for  the  maximum  penetration  depth  in  rock,  concrete,  soil,  ice,  and  marine 
sediments  can  be  found  in  Heuze  (1970),  Zukas  (1990),  etc.  The  maximum  penetration 
depth  is  expressed  as  a  function  of  initial  impact  velocity,  penetrator  cross-sectional  area, 
penetrator  weight,  and  /or  nose  geometry.  However,  in  these  equations  the  only  target 
material  properties  considered  are  the  initial  density  and  unconfined  strength. 

A  large  number  of  semi-analytical  models  where  penetration  is  idealized  as  unifonn 
expansion  of  a  spherical  or  cylindrical  cavity  into  a  semi-infinite  target  have  been  proposed 
(< e.g .  Goodier  (1965),  Hanagud  and  Ross  (1971),  Forrestal  and  co-workers  ( e.g .  Forrestal  et 
al.  (1988))).  These  cavity  expansion  analyses  provide  the  radial  stress  at  the  cavity  surface 
as  a  function  of  the  cavity  velocity  and  acceleration.  This  radial  stress  is  considered  to  be 
the  load  on  the  penetrator  and  is  subsequently  used  in  conjunction  with  equations  of  motion 
to  calculate  the  maximum  depth  of  penetration.  As  emphasized  in  recent  contributions, 
good  agreement  between  model  and  data  can  be  expected  only  if  the  compressibility  and 
rate  sensitivity  of  the  target  material  are  accounted  for  (see  for  example,  Warren  and 
Forrestal  (1998)). 

As  part  of  this  research  effort,  a  new  model  for  calculating  the  resistance  to  penetration  into 
geological  or  geologically  derived  materials  was  proposed.  In  the  analysis  it  is  assumed  that 
the  rate  of  penetration  and  all  flow  fields  are  steady  as  seen  from  the  tip  of  the  penetrator. 
This  hypothesis  is  supported  by  penetration  tests  into  cementitious  materials  (grout, 
concrete)  for  impact  velocities  below  1000  m/s.  The  shear  response  of  the  target  material  is 
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modeled  by  a  modified  Bingham  type  viscoplastic  equation  proposed  by  Cazacu  and 
Cristescu  (2000)  and  Cristescu  et  al.  (2002).  The  rationale  for  adopting  such  a  model  is  that 
it  accounts  for  both  strain  rate  and  compaction  effects  on  yielding,  which  are  key  properties 
of  any  porous/brittle  material.  In  the  low-pressure  regime,  a  non-linear  pressure-volume 
relationship  is  considered.  Based  on  experimental  observations  (e.g.  Schmidt  (2003))  which 
show  that  in  the  high-pressure  regime  a  very  large  increase  in  pressure  is  necessary  to 
produce  even  a  very  small  change  in  density,  the  hypothesis  of  a  “locking  medium”  is 
adopted:  the  density  cannot  exceed  a  critical  value.  The  penetrator  is  wedge-shaped  with  a 
high  length-over-diameter  ratio.  Contact  between  the  projectile  nose  and  target  is 
considered  to  be  of  Coulomb  type  with  constant  friction  coefficient,  whereas  the  frictional 
contact  between  the  remainder  of  the  projectile  and  target  is  considered  to  be  velocity 
dependent.  Resistance  to  penetration  is  calculated  for  different  interface  conditions  between 
the  target  and  wedge.  It  is  shown  that  for  low  to  intermediate  impact  velocities,  accounting 
for  friction  results  in  blunter  optimal  wedge  geometry  for  optimal  penetration  performance. 
Furthennore,  the  higher  the  velocity,  the  greater  is  the  influence  of  nose  geometry  (wedge 
semi-angle)  on  penetration.  Finally,  an  application  of  the  proposed  model  to  a  cementitious 
material  is  presented. 

III.  2.  Constitutive  model  for  the  target 


The  deviatoric  response  of  the  target  material  is  modeled  by  a  non-homogeneous  Bingham 
type  rigid/viscoplastic  equation  proposed  by  Cazacu  and  Cristescu  (2000)  (see  also 
Cristescu  et  al.  (2002)).  The  following  notations  are  used:  (  )  denotes  the  Macauley 


bracket  with  +  fi|)  being  the  positive  part  of  the  real  number^  ,  <r  is  the  Cauchy 


stress  tensor,  S  its  deviator,  /  the  unit  second  order  tensor  ( /y  =  Sy ,  i,  j  =1 . .  .3)  and  p  is  the 
mean  stress  (i.e.  cr  =  £..  +  pS 2  ).  Let  D'  denote  the  deviator  of  the  rate  of  deformation 
tensor  D  defined  as  the  symmetric  part  of  the  spatial  gradient  of  velocity  and  by 
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11  D  = 


2 


D'jD'j  its  second  invariant. 


According  to  the  model  (Cazacu  and  Cristescu  (2000)): 


fo 


D'  = 


if  11  s  ^kfpl 

if  II s  >  k2{p). 


(1) 


The  influence  of  the  degree  of  compaction  on  the  behavior  is  modeled  through  the 
dependence  of  the  yield  limit  in  shear,  k ,  on  the  current  density  p .  Thus,  for  stress  states 

satisfying  II s  <  &(/?)"  the  model  response  is  rigid,  otherwise  the  model  response  is 
viscoplastic.  A  power  law  variation  of  the  yield  limit  with  the  current  density  is  assumed 


k(p)  =  k0  +  p 


f  Y! 

i-^ 

P  ) 


(2) 


where  p0  and  k0  =  k{p() )  are  the  density  and  yield  stress  of  the  undefonned  medium, 
respectively,  while  n  and  ft  are  material  constants.  The  yield  stress  k  is  given  in  stress 
units  and  rj  in  Poise.  The  stress-volume  relationship  adopted  reflects  the  following 
experimental  observations  (see  Schmidt  (2003)): 

•  In  the  low  to  moderate  pressure  regime,  most  cementitious  materials  show  a  highly 
non-linear  mean  stress-volumetric  strain  response,  the  reversible  decrease  in  volume 
being  very  small, 

•  In  the  high-pressure  regime,  a  very  large  increase  in  pressure  is  necessary  in  order 
to  produce  even  a  very  small  change  in  volume. 

Thus,  the  hypothesis  of  a  “locking  medium”  applies,  i.e.  the  density  cannot  exceed  a  critical 
value.  This  critical  density,  called  locking  density,  is  the  density  at  which  no  volume 
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change  occurs  under  hydrostatic  conditions.  It  will  be  denoted  by  p,  and  the  pressure  level 
at  which  this  density  is  first  reached,  called  locking  pressure,  is  denoted  by  pt . 

Hence,  p  =  p*  for  p  >  p,  while  for  p  <  p, ,  the  pressure  vs.  density  relationship  p  =  p  ( p ) 
is  assumed  to  be  of  the  form 


p(p)  =  P* 


(  \ 

P_ 

\P*  j 


(3) 


where  m  is  a  material  constant.  The  relationships  ( l)-(3)  could  be  inverted,  to  give: 


cr  =  < 


p(p)l  + 

p(p)1 


2  rj  + 


Kp± 

4^4 


D'  for  D  *  0, 
for  D'  =  0. 


(4) 


Note  that  the  model  accounts  for  both  rate  dependency  and  compaction  effects  on  yielding, 
which  are  key  properties  of  any  cementitious  or  geologic  materials. 


III.  3.  Computation  of  the  resistance  to  penetration  at  steady-state 
III.  3.1.  Statement  of  the  problem 

Generally,  post-test  observations  of  gun  launch  penetration  experiments  in  geologic  or 
cementitious  materials  (grout,  concrete)  indicate  that  there  are  four  stages  of  the  penetration 
event  (see  Jones  et  al.  (1998)).  The  first  stage  corresponds  to  nose  penetration  when  the 
target  material  is  restrained  from  large  movements;  the  second  stage  corresponds  to 
chipping  and  cratering  at  the  point  of  impact;  the  third  stage  occurs  when  the  projectile  is 
fully  embedded  and  the  pusher  plate  hits  the  target  and  strips  off,  while  the  fourth  stage 
occurs  when  the  full  length  of  the  projectile  is  in  the  penetration  channel  or  tunnel  created 
by  the  penetrator. 

The  analysis  presented  in  this  paper  concerns  the  steady-state  flow  of  the  target  material 
over  a  projectile  fully  embedded  in  the  target  ( i.e .  fourth  stage  of  penetration).  The 
penetrator  is  considered  to  be  rigid  and  wedge-shaped.  Post-test  observations  indicate  that 
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for  impact  velocities  up  to  1000  m/s,  penetration  paths  are  relatively  straight  and  stable 
with  regard  to  the  original  shotline  (see  Jones  et  al.  (1998)).  Thus,  we  can  assume  that  the 
problem  is  axi-symmetric  with  respect  to  the  wedge  centerline.  Let  D  =  D0  U  Dvisco  U  Df 
denote  the  domain  occupied  by  the  penetrator  and  the  target  material  around  it  (see  Figure 
3.1).  The  lateral  extent  of  this  domain  is  F  ■  R  ,  where  R  is  the  wedge  semi-height  and  F  is  a 
number  estimated  from  post-test  observations.  In  the  domains  Do  and  Df  (i.e.  in  front  of  and 
behind  the  penetrator),  it  is  assumed  that  the  target  material  is  in  rigid  body  motion  while  in 
the  domain  DviSCo  the  material  undergoes  viscoplastic  deformation.  We  further  suppose  that 
in  Dvisco  the  flow  lines  are  centered  at  a  certain  pole  O.  Hence,  the  viscoplastic  domain 
Dvisco  is  bounded  by  the  nose  surface,  Tj ,  MP,  and  the  two  surfaces  So  and  Sf.  The  surface  Sf 

R  R 

is  of  radius  OP=rf  = - F ,  while  the  surface  So  is  of  radius  OM  =  r0  =  — - (l  +  F) 

sin  a  sin  a 

(see  Figure  3.1). 


T  N  U 


Figure  3.1:  Schematic  of  the  flow  over  the  wedge. 
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We  assume  that  there  exists  velocity  discontinuities  tangent  to  the  surfaces  S0  and  Sf,  while 
the  components  of  velocity  normal  to  these  surfaces  are  continuous.  Since  across  So 


[4  «=0,  (5) 

with  n  the  outward  normal  unit  vector,  it  follows  that  in  the  viscoplastic  domain,  DviSCo,  the 
only  non-zero  velocity  component  is  the  radial  component, 

Vr  =  V  cos  0  ,  0  <  6*  <a  (6) 

and  consequently  the  only  non-zero  components  of  the  rate  of  deformation  are 

_Vcos6 

r 

and 


Dre  = 


V  sin# 
2  r 


(7) 


The  continuity  equation  reduces  to: 

4~(rP)  =  0 

or 


(8) 


Since  p  =  p0  on  S0,  in  the  viscoplastic  domain, 


P  =  -P„  (9) 

r 

Obviously,  the  compaction  law  (9)  could  also  be  obtained  by  imposing  that  the  velocity  be 
the  same  in  the  domains  Do  and  Df  (debit  compatibility).  In  the  following,  we  estimate  the 
resistance  to  penetration  based  on  energy  considerations.  The  wedge  semi-angle  that 
minimizes  the  resistance  to  penetration  is  then  determined  for  different  interface  conditions 
between  the  wedge  and  the  target. 
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III.  3.2.  Friction  laws 


Along  the  surface  of  the  wedge  significant  friction  arises.  The  mechanics  of  friction  at  high 
sliding  speeds  is  very  complex.  Much  of  the  work  reported  is  at  lower  speeds  or  pressures 
than  those  occurring  during  a  penetration  event.  In  view  of  this,  in  our  analysis  we  assume 
that  a  Coulomb  friction  law  applies  along  Tj ,  the  surface  area  of  contact  between  the  nose 
of  the  rigid  projectile  and  the  target,  i.e. 

T  =  ncrn,  (10) 

where  t  is  the  shear  stress  in  the  deformed  material,  //  is  a  friction  coefficient  taken  as 
constant,  and  crn  is  the  normal  pressure  between  the  nose  and  the  target  material. 

In  plastic  forming  of  metals,  when  one  body  is  fully  plastic  and  the  other  is  rigid,  a  friction 
law  of  the  form  (see  Avitzur  (1968)), 

T  =  jUr<Jy  /  V3  ,  (11) 

is  often  assumed  where  t  is  the  shear  stress  in  the  deformed  material,  pr  is  a  friction 
coefficient  taken  constant  and  <jy  is  the  yield  limit  of  the  material.  In  general,  0  <//,.<  1 , 
where  the  case  //,.  =  0  corresponds  to  “no  friction”,  while  //,,  =  1  corresponds  to  adherence 
of  the  plastic  body  to  the  rigid  die  wall.  Hence,  according  to  the  friction  law  (11) 
irrespective  of  the  normal  pressure  between  the  two  bodies,  the  shear  stress  is  constant.  The 
density  of  the  target  material  behind  the  penetrator  is:  p  =  pf ,  where  p  f -  =  p\r_r  while  the 

state  of  stress  is  such  that =  k(pf)(  i.e.  on  the  yield  surface).  Thus,  it  can  be  assumed 

that  the  contact  between  the  rigid  penetrator  and  the  compacted  material  behind  the 
penetrator  is  described  by  a  law  similar  to  (11),  i.e. 

T  =  prk(pf)  (12) 
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III.  3.3.  Resistance  to  penetration 


The  law  of  conservation  of  energy  is  used  to  compute  the  resistance  to  penetration.  In  the 
domains  Do  and  Df ,  no  deformation  occurs  and  therefore  no  internal  power  dissipation  is 
involved.  In  the  viscoplastic  domain,  the  stress  power  per  unit  volume,  <r :  D  can  be 
expressed  as 

<j  :D  =  ct.D..  =  (477)74,  +  2  k(p)y[u^  +  p-  trD  (13) 

Using  the  work-hardening  law  (2)  in  conjunction  with  the  radial  compaction  law  (9),  we 
obtain: 
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ro  a 

Integrating  over  DviSC0,  we  arrive  at:  W  =  J  0“.DdV=[[(o“.  D)rd6dr  ,  or 
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where  F{a )  =  J  a/3  +  cos2#  d6 . 
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(15) 
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The  surfaces  So  and  Sf  are  surfaces  of  velocity  discontinuity.  The  power  dissipated  at  the 
crossing  of  these  surfaces  is:  JVSg  S/  =  J  t  ■  \y\lA  +  J  t  ■  \y\lA  ,  where  t  is  the  stress  vector 

S0  Sf 

while  [v]  is  the  jump  in  velocity  across  the  respective  surface.  Since  the  magnitude  of  the 
tangential  jump  is  Fsin  9  ,  we  obtain: 


W, 


S0 
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=  J  t  ■  \y\dA  =  J 


<jrg  V  sin  6rQd6  =  - 


rjVz 


■{a  -  sin  a  ■  cos  a )  -  r0  k0 


(hV^)/7,  (a). 


CC  -rr  2 

Ws<  =  •  [v]t£4  =  j-cr/0V  sin  6rfd6  = - ( a  -  sin  a  •  cos  a)  + 
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V  ro  J 


{pS)Fx{a) 


(16) 


where  Ft  (a)=  [  .  =d6 .  Assuming  Coulomb  friction  (see  (10))  along  r, ,  the 

o  v3  +  cos2  6 

surface  area  of  contact  between  the  nose  of  the  rigid  projectile  and  the  target,  the 
dissipation  due  to  friction  is: 


^j,  =  m\(Y  cos  a- age\g=a)di 


r  = 


fiV  cos  a 


—  V  cos  a  +  i)  in 


f  \ 


\rf  J 
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4  cos  a 


3  A 


+  cos'  a 


h{rt,-rr)  +  liff 
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\--t~ 

V  r0  J 


n+ 1 
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(  \ 
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ro 

r0 

-1 

v  P*  ) 

\r f  ) 

(17) 

Assuming  that  the  friction  law  (11)  applies  along  T2 ,  the  surface  area  of  contact  between 
the  penetrator  and  the  compacted  target  material  behind  it  (see  Figure  1),  then 
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(18) 


WTi  =  jur  J  k(pf )  Vds  =  pr  k{p  f )  VL , 

r2 


where  L  is  the  length  of  the  penetrator  shank.  Next,  the  theorem  of  power  expanded  will  be 
used  to  compute  the  resistance  to  penetration.  Consider  the  control  volume 
V  =  D0  U  Dvisco  U  Df  i.e.  in  the  domain  comprised  between  the  planes  x  =  0  and  x  =  OS  and 

y  =  0  and  y  =  F-R,  (see  Figure  1).  At  steady  state,  the  stress  power  theorem  writes: 

J  cr.DdV  +  WSo  +  WSf  +  WTx+Wr2  =  \(<jbV-<jfV)dA. 

D  ■  SV 

VISCO 


(19) 


In  (19),  <jb  denotes  the  magnitude  of  the  stress  vector  acting  along  the  direction  of  the 
velocity  on  the  boundary  of  the  domain  Do  (i.e.  in  front  of  the  penetrator)  while  o  f  is  the 

magnitude  of  the  stress  vector  acting  on  the  boundary  of  Df.  We  may  assume  that 
<7f  =  k{pf),  where  pf  is  the  density  of  the  compacted  target  material  behind  the 

penetrator.  Substituting  (14)  to  (18)  into  (19),  we  obtain  the  expression  for  crb,  the 
resistance  to  penetration,  as: 


V(F  +  l)f? 


(w  +  WSo  +  Ws  +  WTi  +  Wri  +  k(pf)vRF) 


(20) 


III.  4.  Application  to  concrete  (data  after  Schmidt  (2003)) 

To  illustrate  the  predictive  capabilities  of  the  model,  we  will  apply  it  to  concrete.  The 
modeling  effort  was  based  on  the  results  of  a  series  of  quasi-static  compression  tests.  All 
quasi-static  tests  were  conducted  using  a  standard  fully  automated  MTS  testing  machine  at 
a  strain  rate  of  approximately  0.77T0'6/s  and  at  confining  pressures  in  the  range  0-450  MPa. 
(Schmidt  (2003)).  The  density  of  the  material  is  2  g/cm  while  the  uniaxial  compressive 
strength  is  42.2  MPa.  The  parameters  involved  in  the  non-homogeneous  Bingham  model 
(see  eqs.  (l)-(3))  are: 

□  the  viscosity  r) , 
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□  initial  density  p0 , 

□  n  and  /? ,  the  coefficients  involved  in  the  law  of  variation  of  the  yield  stress 
k(p)  (see  (2)), 

□  p„,  p„,  m,  the  parameters  involved  in  the  pressure  density  law  (3). 

For  this  material,  the  locking  pressure  value  is  considered  to  be  p,=  0.5  GPa,  a  value 
which  is  based  on  post-test  data  reported  by  Jones  et  al.  (1998).  The  corresponding  locking 
density  pt  is  estimated  using  the  relationship 


P* 


=  Po 


1  -el 


(21) 


where  s{,  =  0.068  corresponds  to  the  volumetric  deformation  at  failure  under  confinement 
of  450  MPa.  For  sake  of  simplicity,  we  assume  a  quadratic  dependence  of  pressure  on 
density  ( i.e .  we  set  m  =  2  in  (3)).  We  take  the  yield  limit  of  the  undisturbed  medium  k(p0 ) 
to  be  equal  to  the  unconfined  yield  limit,  i.e.  k(p0 )  =  10  MPa  while  the  yield  stress  of  the 
compacted  target  is  taken  to  be  k(pf )  =  FF  k(p0),  FF  being  a  parameter,  which  depends 
on  the  impact  velocity.  Assuming  n  =  2  in  the  law  of  variation  of  k{p )  ,  we  obtain 


P  f  A 

1  _  Po 

{  Pf) 


(22) 


Next,  equation  (20)  is  used  to  calculate  the  resistance  to  penetration.  The  numerical  test 
cases  presented  correspond  to  impact  velocities  in  the  range  of  300m/s  to  1000  m/s.  The 
wedge  dimensions  are  set  at:  R  =  6.35  mm  and  L  =  89  mm.  As  an  example,  Figure  3.2 
shows  (ab/  p, ),  the  resistance  to  penetration  nonnalized  by  the  locking  pressure,  versus  the 
projectile  semi-angle  a  for  an  impact  velocity  V  =  300  m/s.  Calculations  were  done  for 
four  different  values  of  the  coefficient  of  friction  between  the  wedge  tip  and  target  material. 
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The  yield  stress  of  the  compacted  target  k{pf )  was  considered  to  be  2.5  k(p0),  while  the 

coefficient  of  friction  between  the  body  of  the  wedge  and  the  target  was  set  to  p  =  0.6. 
Note  that  irrespective  of  the  contact  conditions  between  the  wedge  and  target,  there  exists  a 
critical  angle  for  which  the  resistance  to  penetration  is  minimum.  Keeping  all  the  other 
parameters  the  same,  the  more  friction  that  is  present,  the  blunter  the  nose  required  to 
achieve  minimum  resistance  to  penetration  (of  course,  assuming  that  no  erosion  takes  place 
and  that  the  nose  does  not  fail).  This  trend  was  observed  for  all  impact  velocities. 


Figure  3.2:  Resistance  to  penetration  crb  nonnalized  by  the  locking  pressure  p*  as  a 

function  of  the  projectile  semi-angle  a  for  different  values  of  the  friction  coefficient 
between  the  target  and  the  wedge  tip;  impact  velocity  V=300  m/s. 
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Optimum  semi-angle  (deg) 


The  plots  of  the  optimal  wedge  angle  vs.  impact  velocity  for  different  values  of  the  friction 
coefficient  between  the  tip  and  target  is  presented  in  Figure  3.3.  Since  compaction  ratios 

were  not  available,  in  the  calculations  it  was  assumed  that  (&(/?y)/ k(p0)  for  an  impact 

velocity  of  1000  m/s  is  double  that  at  300  m/s  impact  velocity.  The  results  indicate  that  at 
higher  impact  velocities  a  sharper  nose  is  required  for  optimum  performance.  Fiowever,  for 
modest  friction  (p  =  0.01),  the  sharpening  is  minimal. 


Impact  velocity  (m/s) 

Figure  3.3:  Optimum  projectile  semi-angle  vs.  impact  velocity  for  various  coefficients  of 
friction  between  the  wedge  tip  and  viscoplastic  target. 
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Figure  3.4  shows  the  resistance  to  penetration  corresponding  to  optimum  wedge  angle 
versus  the  coefficient  //  between  the  wedge  tip  and  target  for  different  impact  velocities. 
Irrespective  of  the  contact  conditions,  an  increase  in  resistance  to  penetration  occurs  with 
increased  impact  velocity. 


0  0.05  0.1  0.15  0.2 

Coefficient  of  friction  tip/target 


Figure  3.4:  Normalized  resistance  to  penetration  corresponding  to  the  optimum  projectile 
semi-angle  vs.  the  coefficient  of  friction  between  the  wedge  tip  and  viscoplastic  target  for 
various  impact  velocities. 
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III. 5  Application  to  concrete  (data  from  Osborn  (1981)) 


The  purpose  of  this  study  is  to  compare  the  steady-state  stresses  on  the  penetrator  obtained 
using  hydrocode  computations  Osborn  (1981)  to  the  steady-state  stresses  obtained  using  the 
proposed  penetration  model.  The  Osborn  (1981)  dry  concrete  model  consists  of  a  hydrostat 
and  a  fixed  yield  surface.  The  hydrostat  defines  a  relationship  between  pressure  and  density 
whereas  the  yield  surface  is  of  Mohr-Coulomb  type.  The  concrete  material  has  an  ambient 
density  p{)  =  2.2  g/cnr  and  an  unconfined  strength  f'c  =  34.47  MPa.  The  hydrostatic 
pressure  (in  GPa)  -volume  relationship  is  given  by: 


'p--\ 

if  0< 

vfo  ) 

lA  J 

14.4 


0.0358  +  7.862 


<  0.0025 


f--1' 

-0.0025 

if  0.0025  < 

\Po  J 

l  Po  ) 
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P--\ 

vA>  J 
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P_ 

vA> 


-1 


<0.3 


(21) 


This  response  curve  is  plotted  in  Figure  3.5. 


78 


CO 

CL 

o 

CD 

i_ 

c/3 

co 

CD 

i _ 

CL 


Figure  3.5:  Hydrostat  for  34.47  MPa  concrete  (after  Osborn  (1981). 


The  yield  surface  is  of  Mohr-Coulomb  type  with  a  saturation  level  set  at  30.67  f'c  .  The 
yield  strength  Y  is  defined  by  the  following  equations: 


'3.3  (p  +  0.1/c73) 
Y  =  \  p  +  2f'c  /  3 

30.67  /; 


if 

f 

if 


-o.i/;/3</7</;/3 
-o.i/;/3<^</;/3 
30.67  /;  <  p 


(22) 
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We  use  (21)  to  determine  the  material  parameters  involved  in  the  pressure-density 
relationship 


p  =  p. 


Indeed,  from  (21)  it  follows  that  the  lock-up  density  is  p,.  =  2.382  g/cm  and  the 
corresponding  lock-up  pressure  is  p,  =  6.3  GPa.  By  least-square  fit,  we  obtain  an  exponent 
m  =12.104.  In  Figure  3.6  are  plotted  together  the  pressure-density  curves  given  by  (21)  and 
(23).  No  data  concerning  the  dynamic  viscosity  nor  the  unconfined  yield  of  the  material 
were  available.  We  assume  the  yield  stress  of  the  undisturbed  medium  k(p0 )  to  be  equal  to 
half  the  unconfined  yield  limit,  i.e.  k(p0 )  =  20  MPa,  set  n  =  2  in  the  power  law  variation  of 


the  yield  stress  with  density  and  that  the  viscosity  is  the  same  as  for  the  concrete  mix 
investigated  by  Schmidt(2003):  rj  =  2- 104  N  s.  The  density  of  the  compacted  region  behind 


T  \  F 

the  projectile  pf  =  p\  _  ,  is  pf  ■  =  —  p0  = - p0 ,  where  2  (F+1)R  is  the  extent  of  the 

r~rf  r.  F 

region  affected  by  the  impact  event.  We  further  assume  F  =  5,  which  results  in  pf  =  2.8584 

g/cm  .  The  dimensions  of  the  penetrator  considered  in  Osborn  (1981)  are:  R  =  1  cm,  L—  14 
cm. 
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Density  (g/cm3) 


Figure  3.6:  Determination  of  the  parameters  involved  in  the  pressure-density  relationship 
(3)  using  Osborn(1981)  data  for  concrete  of  unconfined  strength  of  34.47  MPa 


Figure  3.7  shows  the  resistance  to  penetration  normalized  by  the  locking  pressure  as  a 
function  of  the  projectile  half- angle  for  several  different  striking  velocities  V.  The  yield 
stress  of  the  compacted  target  taken  in  the  simulations  is:k(y»/  )  =  63  k(p0)=  1.26  GPa, 
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the  coefficient  of  friction  between  the  nose  and  the  deformed  target  was  set  to  p  =  0.2 
, while  the  coefficient  of  friction  between  the  shank  and  the  rigid  target  was  pr  =  0.6. 


Figure  3.7:  Resistance  to  penetration  <jh  normalized  by  the  locking  pressure  p*  as  a 
function  of  the  projectile  semi-angle  a  for  different  striking  velocities  V;  k[pf)  = 
63 k(p0),jiu  =  0.2  and  jur=  0.6 


The  plot  of  the  optimal  nose  half-angles  vs.  the  corresponding  impact  velocities  is 
presented  in  Figure  3.8.  Figure  3.9  shows  the  computed  resistance  to  penetration  at 
optimum  angle  for  several  different  impact  velocities.  As  expected,  an  increase  in 
resistance  to  penetration  occurs  with  increased  impact  velocity. 


82 


Optimum  angle  a  (degrees) 


Figure  3.8:  Optimum  projectile  semi-angle  vs.  impact  velocity. 
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Impact  Velocity  (m/s) 


Figure  3.9:  Normalized  resistance  to  penetration  corresponding  to  the  optimum  projectile 
half-angle  as  a  function  of  the  impact  velocity. 


Next,  we  analyze  the  effect  of  interface  conditions.  Figure  3.10  shows  the  optimum  nose 
half-angle  vs.  impact  velocity  curves  corresponding  to  2  different  values  of  the  friction 
coefficients  between  the  nose  and  the  target.  All  the  other  parameters  of  the  model  are  the 
same  as  in  the  previous  simulations.  Figure  3.11  shows  a  comparison  between  the  optimum 
nose  angles  corresponding  to  different  contact  conditions.  Figure  3.12  shows  the  optimum 
projectile  semi-angle  vs.  impact  velocity  for  k{pf)\ &(/0o)  =  63  and  30;  coefficients  of 
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at  optimum  angle  a 


friction://  =  0.2  and  fir=  0.6.  Figure  3.13  shows  the  optimum  resistance  to  penetration 


versus  impact  velocity  curves  corresponding  to  the  yield  stress  ratio  of  63:30. 


Impact  velocity  (m/s) 

Figure  3.10:  Effect  of  friction  between  nose  and  target:  the  lower  the  friction  the  sharper  the 
optimal  nose  angle. 
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at  optimum  angle  a 


Figure  3.11:  Resistance  to  penetration  for  optimum  nose  angle  corresponding  to  different 
friction  conditions  at  the  nose/target  interface. 
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Impact  Velocity  (m/s) 


Figure  3.12:  Optimum  projectile  semi-angle  vs.  impact  velocity  for  k{pf):k(p0)  =  63  :30; 
coefficients  of  friction:  p  =  0.2  and  pr=  0.6. 
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Figure  3.13:  Resistance  to  penetration  crb  normalized  by  the  locking  pressure  p* 
corresponding  to  optimum  nose  semi-angles  as  a  function  of  impact  velocity  for 
k{pr):  k(pQ)=  63  and  30;  coefficients  of  friction:  p  =  0.2  and  pr=  0.6. 
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The  calculated  resistance  to  penetration  values  are  in  good  agreement  with  those  reported 
by  Osborn  (1981). 


III.  5.  Summary  and  Conclusions 

The  purpose  of  this  study  was  the  investigation  of  the  feasibility  of  developing  a  simple 
analytical  model  of  the  planar  steady-state  horizontal  flow  of  a  compressible  viscoplastic 
material  with  locked  hydrostat  between  2  rigid  plates  separated  by  a  distance  2  (F+1)R  and 
symmetric  about  the  penetrator  centerline,  R  being  the  wedge  semi-width.  Resistance  to 
penetration  was  calculated  for  different  interface  conditions.  It  was  shown  that  accounting 
for  friction  alters  the  optimal  wedge  geometry  for  penetration  performance  by  blunting  the 
nose  (wedge  semi-angle).  The  resistance  to  penetration  at  optimal  nose  angle  is  greater  for 
higher  impact  velocity.  Actual  compaction  levels  are  extremely  difficult  to  assess.  Likewise 
is  the  assessment  of  the  friction  that  acts  between  bodies  during  high-velocity  penetration. 
The  analysis  presented  does,  however,  provide  qualitative  insight  into  the  penetration 
process  and  the  role  that  friction  plays  at  low  to  intermediate  impact  velocities.  The 
conclusions  of  this  study  should  be  verified  using  alternative  friction  laws.  This  effort  has 
been  undertaken  in  Cazacu,  Ionescu  and  Perrot  (2004)  where  a  slip-dependent  friction  law 
as  well  as  the  possibility  of  separation  between  the  target  and  penetrator  was  considered. 
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